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Aj  ABSTRACT 

Adsorbed  species  on  a  solid-state  surface  interact  with  the  large  number  of 
substrate  modes,  which  gives  rise  to  thermal  relaxation.  Commonly,  the 
temporal  evolution  of  the  quantum  state  of  the  adsorbate  is  described  by  a 
master  equation  for  the  level  populations  (vibrational  bond,  internal  modes, 
etc).  It  is  pointed  out  that  this  approach  does  not  necessarily  give  a 
correct  account  of  the  coupling  to  the  solid  when  the  effective  level-widths 
become  comparable  to  the  level  separations,  or  larger.  It  is  shown  that  the 
evolution  of  the  populations  does  hot  decouple  anymore  from  the  time 
evolution  of  the  coherences  (off-diagonal  matrix  elements),  which  implies 
that  a  random-phase  approximation  cannot  be  justified,  and  that  the  density 
matrix  of  the  adsorbate  is  not  only  determined  by  the  Golden  Rule  transition 
constants.  Especially  the  line  profiles  turn  out  to  be  very  sensitive  to 
the  coherence -coherence  couplings.  Although  the  coherences  vanish  in 
thermal  equilibrium,  their  time-regression  operator,  and  hence  their  mutual 
couplings  and  their  couplings  to  the  populations,  which  determines  the 
absorption  profile,  does  not.  This  information  is  lost  in  a  master-equation 

treatment  of  relaxation.  , - — - 

fp-  Accession  For  I 


PACS:  71.36,78.90 


I.  INTRODUCTION 


Coupling  of  adsorbed  atoms  or  molecules  on  a  surface  to  the  degrees  of 

freedom  of  the  substrate,  amounts  to  effective  thermal  relaxation  of  the 

internal  molecular  modes,  or  the  vibrational  bond  with  the  solid.  For 

molecules  like  CO  on  a  metal  surface,  the  excitation  of  internal  stretching 

modes  (vibrational  levels)  corresponds  to  a  charge  displacement,  which 

couples  to  the  motion  of  the  electrons  in  the  conduction  band  of  the  metal. 

Interaction  of  the  electron  cloud  of  the  single  molecule  with  the  large 

number  of  electrons  in  the  substrate  then  provides  that  the  small  system 

(the  admolecule)  is  driven  towards  thermal  equilibrium  with  the  heat  bath 

(the  metal). *  This  process  is  usually  regarded  as  electron-hole  pair 

2 

formation  in  the  metal.  Another  example  of  relaxation  of  adsorbates 
pertains  to  the  kinetic  coupling  of  an  adatom  in  a  vibrational 
(electromagnetic)  bond  with  a  crystal  to  the  thermal  motion  of  the  surface 
atoms,  which  support  the  bond.  In  this  fashion  the  kinetic  and  potential 
energy  of  an  atom  in  an  excited  bound  state  on  the  surface,  can  be 

3 

transferred  into  kinetic  and  potential  energy  of  the  crystal  atoms. 

Usually,  this  process  is  viewed  as  a  phonon  exchange  between  the  adatom  and 
the  substrate. ^  Finally,  we  mention  the  irreversible  dipole-coupling  of  dye 
molecules  (coatings)  on  a  dielectric  to  the  surface-plasmon  field. ^  This 
interaction  is  responsible  for  a  dramatic  change  in  the  lifetime  of  the 
molecule,^  which  in  turn  yields  the  celebrated  phenomenon  of  enhanced  Raman 
scattering,  and  surface- induced  (Raman)  resonances. 

Commonly,  the  above-mentioned  mechanisms  for  molecular  relaxation,  due 
to  coupling  with  the  substrate,  are  incorporated  in  the  rate  equations  for 
the  level  populations  n^(t)  (internal  vibrational,  kinetic  or  electronic)  of 
the  adparticle  as 


(1.1) 


k  ni(t)  "  -[(ni(t) 
j 


aij  "  nj^*ji*  +  •  •  •  * 


where  the  transition  constants  a^  are  determined  by  the  Golden  Rule.  Level 
|i>  decays  to  the  other  levels  |j>  at  a  rate  n^a^  (loss  term),  and  gains 
population  from  transitions  of  levels  |j>  to  |  i>  at  a  rate  n^a^.  Balancing 
the  contributions  then  gives  the  master  equation  (1.1),  where  the  ellipsis 

g 

denotes  the  remaining  interactions  with,  for  instance,  a  laser  field.  In  a 

more  thorough  approach,  which  yields  exactly  the  same  result  for  the 

coupling  with  a  heat  bath,  one  starts  with  the  full  Schrodinger  equation  for 

the  density  operator  p(t)  of  the  solid,  the  adparticle  and  the  interaction. 

9  10 

With  standard  reservoir  theory,  projection  techniques  or  a  Langevin 

9 

approach,  one  then  derives  an  equation  of  motion  for  the  reduced  density 
operator  Pg(t)  °f  the  adsorbate,  defined  as 


pQ(t)  *  Trrp(t)  , 


(1.2) 


where  Tr^  indicates  the  partial  trace  over  the  reservoir  states.  Taking  the 
diagonal  part  of  the  equation  then  results  in  the  master  equation  (1.1)  for 
the  populations 


ni(t)  *  < i | pQ( t ) | i>  . 


(1.3) 


If  we  consider  the  off-diagonal  matrix  elements  (the  coherences)  of  p^(t), 
we  obtain 


—  <i|p0(t)|j>  -  {-i(u).  -  -  •|(Ai  +  Aj)}<i|P0(t)| j>  + 


.  ,  (1.4) 


for  i  *  j,  where 


li  »y  a.. 

1  L 
j 


(1.5) 


•  ,  *  r 


V  V  "V..V.V  V  ‘  •*  - •*’. 


equals  the  inverse  lifetime  of  level  |i>,  due  to  transitions  to  every  other 


level  |j>,  and  Kuj.  is  the  energy  of  |i>.  Hence  any  initial  coherence 
< i | Oq ( 0 ) | j >  vanishes  exponentially  on  a  time  scale  (A.  +  A^ )  *  in  the 
evolution  to  thermal  equilibrium  (long-time  limit  t  or  steady  state). 

In  a  random-phase  approximation  one  omits  the  coherences  from  the  beginning 
and  considers  only  the  time  evolution  of  the  populations.  Since  Eq.  (1.4) 
for  the  coherences  does  not  couple  to  the  master  equation  (1.1)  for  the 
populations,  this  approach  is  justified  by  more  elaborated  theories. 

However,  this  is  not  the  whole  story.  First,  the  master  equation  does 
not  embody  the  full  dynamical  behavior  of  the  adparticle.  It  only  provides 
information  on  the  distribution  of  populations,  the  energy  content  of  the 
molecule  and  the  energy  transfer  into  the  solid.  (Eq.  (1.1)  implies  an 
equation  for  energy  conservation,  because  a  level  at  Hw  ,  which  is  populated 
with  n^(t),  corresponds  to  an  energy  n^)W  of  the  molecule.)  Spectral 
information,  like  the  line  profiles  for  the  absorption  of  radiation  and  the 
temporal  correlations  between  the  populations,  which  appear  for  instance  in 
the  response  to  a  pulsed  laser,  are  governed  by  the  time  regression  of  the 
coherences,  rather  than  by  the  master  equation.  Therefore,  in  discarding 
Eq.  (1.4)  one  loses  valuable  information  on  observable  properties  of  the 
system. 

Second,  in  the  derivation  of  the  master  equation  there  is  always  the 


tacit  assumption  that  the  width  of  a  specific  transition  | i>  ■*  |j>,  which 
equals  a^  ,  is  much  smaller  than  the  level  separation  |uu  -  w.  |  .  Then  one 
drops  fast-oscillating  terms  with  the  Bohr  frequencies  4^  *  uk  -  uk  ,  with 
the  argument  that  they  will  approximately  average  out  to  zero.  This  is  a 
correct  procedure  if  |a_|  >>  a„.  In  general,  however,  it  is  not  obvious 
that  this  condition  is  always  satisfied.  To  be  specific,  for  electron-hole 


-1 


acquire  values  of  the  order  of  a 


pair  formation,  both  a,,^  and  A,. 

ij  ij 

picosecond,  and  for  CO  on  copper,  for  instance,  the  inverse  lifetime  of  the 
stretching  modes  exceeds  the  level  separations.  In  the  case  of  adsorption 
of  atoms  on  a  crystal,  the  timescales  for  relaxation  typically  range  from 
nanoseconds  to  picoseconds.  For  weakly- bounded  atoms  (shallow  potentials) 
the  levels  become  closely  spaced,  and  their  separations  can  easily  assume 
the  order  of  magnitude  of  the  damping  constants. 

In  this  paper  we  shall  retain  the  couplings,  which  are  due  to  the 
overlap  of  levels,  and  point  out  their  significance.  Especially  the 
absorption  profile  for  weak  radiation  appears  to  be  sensitive  to  the 
coherence-coherence  coupling,  as  will  be  illustrated  by  an  example. 


II.  RELAXATION  THEORY 


In  order  to  display  clearly  the  various  approximations  and  to  set  the 

notation,  we  summarize  the  basic  elements  of  relaxation  theory  in  this 
12 

section.  With  Ha  the  Hamiltonian  of  the  adsorbate,  which  includes  the 

binding  energy  to  the  surface  (attractive  force  of  electromagnetic  origin), 

and  Hf  the  substrate  Hamiltonian,  the  equation  of  motion  for  the  density 

operator  p(t)  of  the  active  system  reads 

irt^r  p(t)  -  [H  +  H  +  SR,  p( t ) ) ,  p(t)+  *  p(t),  Trp(t)  =  1  .  (2.1) 

at  a  r 


Here  the  interaction  between  the  molecule  and  the  substrate,  which  gives 

rise  to  the  relaxation,  is  divided  as  SR,  where  S  *  S+  is  a  molecular 

operator  (S  =  systems),  and  R  =  R^  (reservoir)  is  an  operator  in  the 

subspace  of  the  solid.  For  instance,  for  coupling  by  single-phonon 

transitions,  the  operator  S  equals  the  derivative  of  the  binding  potential 

well,  perpendicular  to  the  surface,  and  R  is  the  amplitude  operator  of  an 

1 3 

atom  in  the  crystal.  With  a  Taylor  expansion  it  is  always  feasible  to 


factorize  the  interaction  as  but  in  order  to  avoid  many  obfuscating 

subscripts,  we  only  retain  a  single  term.  In  the  case  of  phonon  coupling, 
the  subscript  i  indicates  the  number  of  phonons  which  are  involved  in  a  one 
step  transition.  By  writing  SR  for  the  interaction,  we  only  keep  track  of 
single-phonon  transition.  This  is  already  very  accurate  if  the  transition 
frequencies  are  smaller  than  the  Deybe  frequency  of  the  crystal,  since 
then  every  pair  of  levels  is  resonantly  couplied  by  a  single-phonon 
interaction. 

The  standard  integral  of  Eq.  (2.1),  which  is  suitable  for  the 
development  of  relaxation  theory,  is  most  conveniently  derived  with  the  aid 
of  the  interaction  picture.  With  the  Liouvillians 

L  •  =  h_1[H  ,•]  ,  L  •  =  h_1[H  ,•)  ,  (2.2) 

a  a  r  r 

we  define  the  transformed  density  operator  as 

o(t)  =  exp( i(L  +  L  )t)  p(t)  ,  (2.3) 

3  r 

which  obeys  the  equation  of  motion 

ih4ro( t )  =  [ S( t )R( t ) , o( t ) ]  .  (2.4) 


The  free  evolution  of  S  and  R  is  displayed  in  their  time  dependence, 
according  to 

S(t)  -  expfiL^OS  ,  R(t)  =  expfiL^tlR  .  (2.5) 

Iterating  Eq .  (2.4)  twice  and  differentiating  the  result  with  respect  to 
tine  then  yields  the  integral 

i^-o(t)  =  h'1[S(t)R(t),o(0)]  { 


Jo 


[S(t )R( t ) , [S(t,)R(t,)fo(t,)J] 


(2.6) 


As  the  initial  value  we  can  choose  arbitrarily  o(0)  =  oft(0)p  ,  with  o_(0)  = 

u  eq  U 

Tr  o(0)  the  reduced  density  operator  of  the  adsorbate  and  p  the  thermal- 
r  eq 

equilibrium  state  of  the  solid.  By  definition,  pg^  commutes  with  Hf,  and  we 
shall  assume  that 


Tr  Rp  =  0  .  (2.7) 

r  eq 

This  means  that  the  interaction  does  not  give  rise  to  a  net  energy  content 

of  the  molecule,  which  is  exact  in  most  cases.  Otherwise  the  restriction 

(2.7)  can  be  enforced  by  a  proper  transformation  of  p(t),  which  shifts  the 

interaction  over  its  average.  As  a  second  step  we  write  o(t')  =  o.(t')p  + 

0  eq 

correction,  which  can  always  be  done.  After  taking  the  trace  in  Eq.  (2.6) 
over  the  reservo: r  states,  we  then  obtain 

dt°0( t )  =  *  2tt  {  dT{G(T)S(t)S(t-T)oQ(t-T) 

+  G(t)  Og(t-T)S(t-T)S(t)  *  G(t)  S(t)oQ(t-T)S(t-T) 

-  G(T)S(t-t)oQ(t-t)S(t))  +  correction  .  (2.8) 

For  later  purposes  we  do  not  impose  the  condition  oj  =  Oq,  which  would  allow 
a  contraction  of  term^ .  The  interaction  with  the  substrate  is  incorporated 
in  the  complex-valued  reservoir  correlation  function 

G( t )  =  2trK'2Tr  R(t)Rp  ,  (2.9) 

r  eq 

which  can  be  evaluated  immediately  for  any  prescribed  R  and  H^.  Among  other 
features,  this  G(x)  encompasses  the  temperature  dependence  of  the  coupling. 
Equation  (2.8)  is  the  usual  starting  point  for  reservoir  theory. 


A  particularly  transparent  representation  of  Eq.  (2.8)  emerges  if  we 
adopt  a  Liouville-operator  notation  for  the  couplings  in  the  integrand.  To 
this  end,  we  substitute  S(t)  =  exp(iLat)S,  transform  back  to  the  p-picture, 
and  define 


V  *  ts.o 


Then  Eq.  (2.8)  assumes  the  form 


4»0(t)  ■  Vo(t)  -  2*S  aT<Sp„(t-,» 

A  -iL  t 

-G(t)  e  3  (pQ(t-T)S)}  +  correction  . 


(2.10) 


(2.11) 


Next  we  introduce  the  correlation  operator 

L  (t)  •  =  S-G(t)  -  G(t)**S  ,  (2.12) 

c 

which  enables  to  rewrite  Eq.  (2.11)  compactly  as 

l5rp0(t)  -  W0  •  k  LS  6  S  Lc(T)p0(t-’)  +  C0rreCtl0n  ■ 

(2.13) 


As  long  as  the  correction  terms  are  taken  into  consideration,  and  provided 
that  condition  (2.7)  holds,  this  is  still  an  exact  integral  of  Eq.  (2.1). 
Equation  (2.13)  clearly  reveals  the  significance  of  L  (t).  It  accounts  for 
the  memory  in  the  interaction. 


III.  RESERVOIR  CONDITION 

A  solid  has  a  broad,  continuous  spectrum  of  modes,  which  are  coupled  to 
the  molecule  by  the  SR  interaction.  This  property  assures  that  the 
correlation  fu.  ction  G(t)  from  Eq.  (2.9)  decays  to  zero  very  fast  for 


increasing  t.  Typically,  the  time  width  of  G( t )  is  of  the  order  of  the 
inverse  cut-off  frequency  of  the  mode  distribution,  which  is  the  Debye 
frequency  for  a  harmonic  crystal.  Hence  this  time  is  much  shorter  than  any 
tine  scale  on  which  Og(t)  varies,  due  to  its  coupling  with  the  same 
continuum.  If  we  indicate  the  correlation  time  of  the  reservoir  by  t  ,  and 
the  damping  constant  for  the  transition  |  i>  -*■  |j>  by  a^j*  as  Eq.  ( l  - 1 ) » 
then  the  reservoir  condition  reads 


a .  .  t  <<  1 
ij  c 


(3.1) 


for  every  i,j.  For  rare-gas  atoms  on  a  crystal  like  KC1,  typical  values  of 
the  product  ajjTc  ran8e  from  0.1  to  0.01.  By  definition,  the  value  of  the 

integral  in  Eq.  (2.13)  has  an  order  of  magnitude  of  and  it  is  easy  to 

12  2 
estimate  that  the  correction  terms  are  of  the  order  of  a. ,  t  .  Therefore 

lj  c 

they  can  be  discarded  in  comparison  with  a_.  In  the  interaction  picture, 
o^t-t)  evolves  on  a  time  scale  a^  *,  which  implies  that  we  can  replace 
on(t-t)  by  <Jn(t)  in  Eq.  (2.8).  In  the  p-picture  the  density  operator 


oscillates  with  the  Bohr  frequencies,  and  hence  the  corresponding 


approximation  in  Eq.  (2.13)  is 


.  14 


exp( iL  (t-r))pn(t-T)  -  exp( iL  t)p.(t)  . 

a  u  “  a  u 


(3.2) 


Then  we  can  take  Pq(0  outside  the  integral,  and  again  with  condition  (3.1) 
we  can  replace  the  upper  integration-limit  by  t  =  ®.  Combining  everything 
then  leads  to  the  concise  form  of  the  equation  of  motion 


The 


P0 


(t)  =  (L  - 
a 


Liouville  operator 


ir)pQ(t)  . 

r  is  defined  by 


-iL  t  iL  t 

r=r— L  die  aL  (t)e  3 

2"  s  J0  c 


(3.3) 


I 


(3.4) 


which  accounts  for  the  relaxation  of  the  molecule,  due  to  the  coupling  with 


the  modes  of  the  substrate.  Recall  that  the  only  assumption  in  deriving  Eq 
(3.3)  is  the  condition  (3.1),  which  is  fairly  justified  for  adsorbates. 
Equations  (3.3)  and  (3.4)  constitute  the  full  effect  of  the  coupling  to  the 
substrate,  including  the  time  evolution  of  the  coherences,  coupling  between 
coherences,  level  shifts,  etc.  The  advantage  of  the  Liouville  approach,  in 
contrast  to  a  master -equat ion  treatment,  which  relates  matrix  elements,  is 
that  the  dynamics  of  the  interaction  is  completely  disentangled  from  the 
structure  of  the  equation  of  motion,  Eq .  (3.3). 


IV.  EVALUATION  OF  T 


Equation  (3.4)  defines  the  relaxation  operator  F  in  terms  of  the 


Liouvillians  L  and  L  (t).  If  we  remember  that 
a  c 


exp(?iL  t)-  =  exp(tiH  i/h)-  exp(+iH  i/ft)  , 

3  3  3 


(4.1) 


and  then  insert  the  definition  (2.12)  of  L^x)  into  Eq .  (3.4),  we  find  that 


F  can  alternatively  be  represented  by 


-  L  (Q-  -  -Q  )  , 


(4.2) 


which  involves  the  Hi lbert -space  operator 


■fcC 


•  iL  x 


dxG( t  )e 


(4.3) 


Notice  that  Eq .  (4.2)  still  contains  four  terms,  because  L  is  a  commutator 

s 


This  also  immediately  implies  the  properties 


( To ) +  =  ro+  ,  Tr  (To)  =  0  , 


(4.4) 


for  any  a.  These  relations  are  necessary  for  the  conservation  of 
Herroiticity  and  trace  in  the  time  evolution  of  pg(t). 

Since  appears  in  an  exponent  in  the  definition  of  Q,  an  expansion  in 
matrix  elements  is  most  convenient  on  a  basis  where  is  diagonal,  e.g., 

the  adsorbate  states  |i>.  For  simplicity  we  shall  suppose  that  the  states 
|i>  are  non-degenerate.  This  is  no  restriction  at  all,  but  it  avoids 
cumbersome  notations.  In  terms  of  the  projectors 


P.  =  |i><i|  , 

the  exponential  in  Eq.  (4.3)  can  be  expanded  as 

p(-iL  t)-  =  )  exp(-iA. ,t)P.  ■  P. 

F  a  /.  ij  1  J 


(4.5) 


ext 


(4.6) 


ij 


Then  Q  assumes  the  simple  form 


Q  =  ~  j  g.  .P.SP.  , 

2  L  J1  1  J 


(4.7) 


ij 


where  the  reservoir  parameters  g  are  defined  by 


8ji  =  ^dT  eXp('i4ijT)G(T)  ’ 


(4.8) 


which  is  essentially  the  Laplace  transform  of  G(t).  Substituting  Q  and  its 
Hermitian  conjugate  into  Eq.  (4.2)  gives 

r°  »  J-  l  S..{gjiS|i><j|o  -  g.*  SojiXj| 

ij 


-gjJiXjloS  +  g.j  °|  i><j  I S }  , 


(4.9) 


which  defines  the  action  of  T  or  an  arbitrary  Liouville  vector  o.  Here 


< i | S | j >  is  a  matrix  element  of  the  molecular  part  of  the  interaction. 


' '"Jl 


"J  ,mJ> 


■  F  *"  T— "W  -  W  -  •  ”  V  -  V  -  V  * 


Equation  (4.9)  is  a  hybrid  representation  of  r,  since  it  contains 


matrix  elements  ,  as  well  as  the  operator  S  itself.  With  the  closure 


relation 


h 


(4.10) 


we  can  cast  Eq.  (4.9)  in  the  form 

Fo*2  1  {ckiijlk><jl°  +  Ck!ij°|j><k)} 


ijk 


2  \  <ciaijli><j|°!k><tl  +  ckiijlt><kl°|j><il>  • 

ijkt 


(4.11) 


where  we  introduced  the  coefficients 


C  as  S  S  S  • 

ktmn  kt  mnRnm 


(4.12) 


Hence  the  relaxation  operator  is  determined  by  the  reservoir  parameters  g 


iJ 


and  the  matrix  elements  S..  between  the  wave  functions  I i>  and  |j>  of  the 

XJ  I  I J 


adsorbate.  Both  g. .  and  S. .  can  be  evaluated  directly,  once  the  kind  of 
ij  ij 


substrate  and  the  molecular  wave  functions  (internal  modes),  or  the  binding 
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potential  (vibrational  coupling  to  a  crystal),  are  prescribed. 

The  relaxation  gives  rise  to  an  effective  width  of  the  levels  (their 
inverse  lifetimes).  If  we  assume  that  the  overlap  between  different 
resonances  is  negligible,  then  we  only  have  to  retain  couplings  between 
pairs  of  levels  in  Eq.  (4.11),  rather  than  between  four  levels 


simultaneously.  If  we  further  neglect  the  imaginary  part  of  g„  (the  level 
shifts),  then  Eq.  (4.11)  attains  the  familiar  form^ 


To  =  ^  ^  a^j{P^o  +  oP.  -  2Pj  <i|o|i>)  , 


(4.13) 


ij 


,  r,  v. 


-  .V 


12 


l_V  JV.J*. 


•  ***•'•'  *  •  '  *.  .  %  4 , '  %  ‘ .  i  \  '»  .%  .*»  -**  , 


& 


y: 

s.‘ 


'£n 


,v , 


•>i 


& 


tjCm 


in  terms  of  the  rate  constants  a_  for  the  transitions  |  i>  •*  |j> 


a.  .  =  Re  c.  =  Is,  .1  Reg,  . 
ij  kttk  1  kfc'  6kt 


(4.14) 


It  is  easy  to  check  that  the  simplified  form  (4.12)  leads  to  the  master 
equation  (1.1). 


V.  EQUATION  OF  MOTION 

Numerical  solution  of  the  equation  of  motion  (3.3)  requires  an 
expansion  in  matrix  elements.  With  Eq.  (4.11)  for  the  relaxation,  and 


Lap0  =  \  wi(li><ilpo  ”  poli><il)  * 


(5.1) 


for  the  free  evolution,  we  obtain 


^<m|p0|n>  =  -iAmn<m|p0|n>  -  ±-  ^..Xj  |  PQ  |  n>  +  cn.*.<n»|  pQ  |  j» 

ij 

+  2  I{cinmj<jlpoli>  ^  ciInj<ilpolj>}  ‘ 

ij 

If  we  subsequently  set  m  =  n  and  use  p^  =  p^,  we  find 

^<m|p0|m>  =  -  I  ^{cmiij<j|p0|m>  -  c  .^j  |  pQ|  i» 
ij 

+  Hermitian  conjugate 


(5.2) 


(5.3) 


for  the  time  derivative  of  the  population  <m|pg|m>.  This  result  is  by  no 
means  equivalent  to  the  master  equation  (1.1).  Even  if  we  were  to  neglect 
the  imaginary  part  of  the  coefficients,  then  the  time  evolution  of  the 
populations  would  still  couple  to  the  real  part  of  the  coherences  < j | p^ J m> 
(j  *  m),  due  to  overlapping  resonances.  Conversely,  the  time  evolution  of 
the  coherences,  as  it  follows  from  Eq.  (5.2)  with  m  *  n,  couples  to  the 


populations  and  to  every  other  coherence,  in  contrast  to  the  simple 
exponential  decay  (Eq.  (1.4))  in  the  limit  of  non-overlapping  resonances. 

In  conclusion,  if  the  widths  of  transitions  become  comparable  with  the  level 
separations,  it  is  inevitable  to  take  the  coherences  into  consideration, 
which  in  turn  prohibits  the  derivation  of  a  genuine  master  equation.  Stated 
otherwise,  in  this  case  a  random-phase  approximation  cannot  be  justified. 


VI.  STEADY  STATE 

Of  particular  importance  is  the  long-time  solution  p^  =  Pg(t  ■*  ®). 

This  steady-state  density  operator  obeys  dp^/dt  =  0,  and  its  matrix  elements 
follow  from  Eqs.  (5.2)  and  (5.3)  with  the  left-hand  sides  replaced  by  zero. 
The  relaxation  of  the  reduced  density  operator  to  a  steady  state  is  a  result 
of  the  coupling  to  the  substrate.  Apart  from  extravagent  situations,  this 
Oq  is  unique.  Therefore  we  can  solve  Eq.  (5.3)  by  trial  and  error.  If  we 
assume  that  the  coherences  vanish  in  the  steady  state,  e.g. 


<  i  I  I  j  >  *  6..<i|pn|i>  =  <5..n. 

■  0,J  i]  1  O'  lj  l 


then  Eq.  (5.3)  reduces  to 


)(n.a..  -  n.a..)  =0  ,  )  n.  =  1  , 

l  i  ij  j  j  1  L  1 


(6.1) 


(6.2) 


where  the  rate  constants  a„  are  defined  in  Eq.  (4.14).  This  equation  is 
exactly  the  same  as  Eq.  (1.1)  for  t  •*  ®,  so  with  the  left-hand  side  set 
equal  to  zero.  Hence  the  long-time  solution  of  the  density  matrix  is 
diagonal,  and  determined  by  the  Golden  Rule  transition-constants  a„.  Note 
that  the  imaginary  parts  of  c^^  do  not  contribute.  In  terms  of  the 
solution  n^  of  Eq.  (6.2)  and  the  projectors  P^,  we  can  write  Pq  as 


50  *  l  Vk  • 


(6.3) 


-  Jb.j*  \  ■>.1kkaa  kauWAVVWA’A  ^ 


which  will  be  used  in  the  next  section. 

Only  in  the  transient  region  the  coherences  affect  the  details  of  the 
tine  evolution  of  the  populations  n^(t).  Therefore  it  might  appear  that  the 
importance  of  the  coupling  to  the  coherences,  and  in  particular  the  coupling 
between  different  coherences,  is  marginal.  Such  is  however  not  the  case. 
Even  in  the  steady  state,  where  the  coherences  disappear  identically,  their 
time  regression  has  a  great  significance  for  the  calculation  of  observable 
quantities.  It  should  be  emphasized  that,  for  instance,  a  level  population 
is  not  directly  accessible  to  experimental  observation  in  general. 

VII.  ABSORPTION  PROFILE 

A  common  method  to  observe  resonances  of  adsorbates  is  by  irradiation 
of  the  surface  with  a  low-power  monochromatic  laser,  and  measuring  the 
absorption  as  a  function  of  the  photon  energy  (frequency).  Probing  the 
system  with  a  weak  radiation  field  has  the  advantage  that  it  does  not 
disturb  the  molecule  (excitation  of  internal  modes),  nor  desorb  the  layer  or 
heat  the  substrate.  Furthermore,  the  intensity  of  the  radiation  (scattered 
or  transmitted)  can  be  detected  with  a  high  accuracy,  partly  due  to  the  fact 
that  only  a  relative  measurment  is  required,  which  relates  the  absorption  at 
a  particular  frequency  to  a  calibrated  off-resonance  background  level.  Care 
should  be  exercised,  however,  in  the  interpretation  of  the  spectrum.  It  is 
not  the  free  molecule  which  absorbs  the  radiation,  but  the  joint  system  of 
molecule,  substrate  and  interaction.  Conversely,  this  feature  provides  an 
interesting  tool  to  investigate  the  molecule-surface  interaction  (binding 
potential,  charge-exchange  mechanism),  or  properties  of  the  substrate 
(dispersion  relations  of  phonons  or  polaritons).  Especially  for  a 
transparent  crystal,  like  for  instance  IR  light  on  KC1,  information  about 
the  crystal  can  be  obtained  by  spectroscopic  methods  in  this  way,  which  is 
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i 

not  feasible  without  the  coating.  The  adbond-mediated  absorption  profile  j 

will  reflect  the  details  of  the  coupling  mechanism  of  the  molecule  to  the 

L 

substrate,  and  the  properties  of  the  solid  itself.  We  remark  that  even  if  j 

both  the  adsorbate  and  the  substrate  are  transparent  for  a  particular 

frequency,  the  bond  can  be  optically  active.  This  is  for  instance  the  case 

for  rare-gas  atoms  on  an  alkali-chloride  crystal.  Then  the  absorption 

profile  will  reveal  the  properties  of  the  wave  function  of  the  vibrational 

state,  and  as  we  show  below,  in  great  detail. 

If  we  indicate  by  p  the  dipole-moment  operator  of  the  molecule  or  the 
bond  (or  both),  and  by  e  the  polarization  of  the  probe  beam,  then  the 
absorption  profile  is  determined  by  the  time  regression  of  the  operator 

d  =  p- e  .  (7.1) 

With  d(t)  the  Heisenberg  representation  of  d,  the  number  of  absorbed  photons 
per  unit  of  time  from  the  incident  field  with  frequency  u  equals 

IU)  =  I  (cnh2c)*1Re  fdt  e1U)t  Tr  p[d(t)+,d]  ,  (7.2) 

P  0  J0 

where  1^  is  the  laser  power  (energy  per  unit  of  time  through  a  unit  area, 
perpendicular  to  the  direction  of  propagation).  Multiplication  of  I(uj)  by 
gives  the  absorbed  energy  per  unit  of  time.  Here,  p  is  the  thermal- 
equilibrium  state  of  the  entire  system,  and  the  time  regression  of  d(t)  is 
governed  by  the  Hamiltonian  +  SR  from  Eq.  (2.1). 

After  a  transformation  to  the  Schrodinger  picture,  we  can  eliminate  the 
reservoir  degrees  of  freedom  from  the  integrand  in  Eq.  (7.2)  along  exactly 
the  same  lines  as  we  derived  Eq .  (3.3)  for  the  reduced  density  operator 


pg(t).  We  obtain 


-i(L  -ir)tf 

Trp[d(t)+,d]  »  Tr  d+e  [d,p0J  * 


where  Tra  runs  over  the  adsorbate  states  only.  Equation  (7.3)  can  be 

regarded  as  Liouville-operator  representation  of  the  quantum-regression 
1 8 

theorem.  Then  we  insert  Eq.  (7.3)  into  Eq.  (7.2)  and  perform  the  time 
integration,  which  yields  the  representation 

«“>  -  tp<c/c>'W/  -  rrTTr  ld’5o]  (7-‘ 


in  terms  of  an  operator  inversion.  The  upper  limit  t  *  °»  does  not 
contribute,  due  to  the  identity 


-  i(L  -iDt 

lim  e  3  [d,pQ]  *  p0  Tra^d,Po'  *  °  * 

t-~°  J 


which  in  turn  follows  from  the  facts  that  pQ  is  the  solution  of 


(7.5) 


(L  -  ir)pn  *  0 
a  u 


(7.6) 


and  that  the  trace  is  conserved  in  a  time  regression  with  exp( - i( L^- i r )t ) . 
We  recall  that  the  commutator  ( d , Pq ]  gives  rise  to  the  terms  dp^  and  p^d , 
which  correspond  to  stimulated  absorption  and  emission  in  the  laser  mode 
respectively.  The  profile  I(ai)  is  the  balance  between  these  two  processes 
Since  we  suppress  the  degeneracies  of  the  levels,  we  can  expand  the 
operator  d  on  the  adsorbated  states  as 


d  -  l  dijltxJl  • 


(7.7) 


With  the  representation  (6.3)  of  the  steady-state  solution,  we  find 


[d*P0l  =  [(nj  ‘  ni)dijli><j|  • 


(7.8) 


Substitution  into  Eq.  (7.4),  expanding  d^  as  in  Eq.  (7.7),  and  taking  the 


trace  then  gives 


■<»>  *  V'o*2^'1  l  (Sj  •  fii)Re  diA*<k|<*  -  rr  irli><jl)lt>  • 

ijkl  a 


(7.9) 


The  matrix  elements  of  w  -  +  if  can  be  read  off  immediately  from  Eq. 

(5.2),  and  after  inversion  of  this  matrix,  Eq.  (7.9)  yields  the  absorption 

profile.  Hence  for  a  system  of  N  levels,  the  evaluation  of  I(oj)  requires 

2  2 

only  the  inversion  of  an  N  *  N  matrix. 

From  the  explicit  representation  (7.9)  it  follows  that  I(w)  is 
determined  by  the  time  regression  of  the  coherence  | i>< j | ,  (i  *  j).  In  the 
tine  domain,  like  in  Eq .  (7.1),  we  propagate  |i><j|  with  exp( - i ( L  - iT )t ) , 
and  then  take  the  coherence  <k|...|l>  of  the  result  at  time  t.  A  Fourier 
transform  then  gives  the  spectral  profile.  This  elucidates  the  importance 
of  the  coherences  in  the  observation  of  adsorbates  with  spectroscopic 
methods.  It  is  their  time  regression  which  determines  the  spectral 
distribution,  whereas  the  populations  only  appear  as  an  overall  factor  n^  - 
n  .  Furthermore,  the  spectral  resolution  involves  the  complete  operator  L ^ 

-  if,  which  represents  the  free  evolution  of  the  molecule,  the  damping,  the 
level  shifts,  and  all  the  couplings  between  coherences  and  populations. 

In  the  secular  approximation,  Eq.  (4.13),  where  only  couplings  between 
pairs  of  levels  are  taken  into  account,  only  the  terms  with  i  *  k,  j  *  t 
survive  in  Eq.  (7.9).  We  notice  that  for  (i,j)  *  (k,fc)  an  overall  factor 

*  i  1 2 

ar*ses»  rather  than  |d.^|  .  This  product  of  matrix  elements  carries 

information  on  the  relative  phases  of  the  dipole-moment  matrix  elements  d. . , 

2 

which  is  not  the  case  in  the  secular  limit,  where  only  |d^|  appears. 
Therefore,  for  transitions  which  have  a  sufficient  overlap,  it  should  be 
feasible  to  extract  their  relative  phase  from  an  absorption  measurement. 


VIII.  TWO-LEVEL  SYSTEM 


In  order  to  display  the  significance  of  the  various  notions,  we  utilize 
the  two-level  system.  For  the  case  of  a  vibrational  coupling  of  an  atom  to 
a  crystal,  this  situation  can  occur  if  the  potential  well  is  shallow  enough, 
so  that  it  supports  only  two  discrete  levels.  The  steady-state  populations 
are  readily  found  from  Eq.  (6.2).  We  obtain 
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(8.1) 


Next  we  neglect  the  self-coupling  of  a  level,  which  means  that  we  assume 
S 1 1  =  S^2  =  0.  The  equation  for  the  coherence  p^^(t)  follows  from  Eq. 
(5.2).  Explicitly 


dt  p21(t)  =  'lA21p21(t)  '  2(c21 12  +  C1221)p21(t) 


+  2(c2121  +  Cl2l2)p12(t)  ’ 


(8.2) 


and  the  equation  for  p^t)  =  P21(t)  f°H°ws  after  a  complex  conjugation. 


We  shall  write  =  "q  ^or  t*ie  resonance  frequency  of  the  adsorbate,  before 


coupling  to  the  reservoir,  and  we  introduce  the  notations 


H  =  2 1 S21 1 2(g21  +  812}  ’  (8‘3) 

S21  "  |S2i  |exp(i<>)  ,  (8.4) 


where  the  latter  defines  the  phases  $  of  Then  Eq.  (8.2)  attains  the 

form 

^•P21(t)  =  -( iojQ  +  n)P21(t)  +  n  exp(2i0)p12(t),  (8.5) 

which  shows  immediately  that  =  p^  =  0,  in  agreement  with  Eq.  (6.1). 
The  second  term  on  the  right-hand  side  represents  a  coupling  between  the 


coherences  p^^Ct)  and  p^vt).  This  term  would  be  discarded  in  a  secular 
approximation. 

Equation  (8.5)  together  with  its  complex  conjugate  constitutes  a  closed 
set,  which  has  the  characteristic  frequencies 

U)+  *  t{<u0(oj0  +  i(n  -  n))  -  4^21  +  al2^2^  '  2^a21  +  a12^’  ^8-6^ 


where  we  used 


n  ♦  n  -  »2l  *  *l2  • 


(8.7) 


In  the  secular  limit  tu^  >>  a^  +  a^>  t^iey  reduce  to  w+  -  i  ujq  -  iiCa^^  + 
a^),  and  hence  the  resonance  is  situated  at  uj^  with  a  width  equal  to  i(a2^ 

+  al2^'  Conversely,  for  a^^  +  a^  >>  we  find  w+  -  0,  w_  -  -i(a2j  +  a^). 
which  are  both  situated  at  id  =  0. 

The  full  absorption  profile  I(id),  Eq.  (7.9),  is  determined  by  the 
Liouville  operator  u>  -  L  +  iT.  From  Eq.  (5.2)  we  can  easily  construct  a 

3 

matrix  representation  for  -i(La  *  if),  which  in  turn  gives 


ui  +  ia_ 


w  +  ia. 


ui  -  L  +  i  r 
a 


w  *  Wq  +  in 

- inexp( -2 i$) 


-in  exp(2i4>) 


ui  +  u)q  +  in 


(8.8) 


on  the  basis  |2><2|,  |l><l),  |2><l|,  |l><2|.  Inversion  of  this  matrix  and 
substitution  into  Eq.  (7.9)  then  readily  yields 


p  a^^  +  a^  11  (w  +  u)q  +  in*)(w  -  <*>q  +  in)  +  nn 


where  we  introduced  the  Einstein  B-coeff icient 


pertaining  to  stimulated  transitions  between  |l>  and  |2>.  Notice  that  the 


phase  0  of  disappears  from  I(id),  as  could  be  anticipated.  Only  relative 
phase  differences  between  wave  functions  of  pairs  of  levels  might  have  a 
significance,  and  hence  the  phase  difference  $,  related  to  the  wave 
functions  of  a  two-state  system,  should  vanish  in  observable  quantities. 

The  resonances  of  I ( o>)  are  situated  at  the  real  parts  of  the  zero's  u>+  of 
the  denominator.  Since  Re  id+  =  -Rein  there  is  only  a  single  peak  in  I(id) 
for  positive  frequencies  id. 

In  the  secular  limit  we  can  omit  the  term  nn*  in  the  denominator  (with 
2 

respect  to  id^  ),  and  then  the  profile  reduces  to 


IU)  =  BI 

s  pa 


21  a12  1  „ 

- - - Re  — 

21  +  a12  "  w 


“0  +  in 


(8.11) 


which  is  a  Lorentzian  around  Wq  +  Iran,  with  a  half  width  at  half  maximum 

equal  to  Re  n  =  i C a ^ ^  +  a^)-  Figs.  *  ar)d  ^  we  ^ave  Plotted  I(cd)  and 

I(~)  for  Im  n  =  0  and  for  different  values  of  Re  n-  Note  that  I(id)  and 
s 

I ( — )  assume  the  same  value  at  the  transition  frequency  id-,  if  Im  n  =  0. 
s  u 

Increasing  Re  n  shifts  the  peak  in  I(id)  towards  lower  frequencies,  without  a 
significant  broadening.  On  the  other  hand,  an  increas  of  Ren  in  I(id)s 
broadens  the  line,  without  shifting  it.  The  qualitative  different  effects 
on  I(id),  as  compared  to  the  approximation  I(w)  ,  arise  purely  due  to  the 
inclusion  of  the  coherence-coherence  coupling.  The  shift  of  the  line  in 
I(_)  should  not  be  confused  with  the  shift  Im  n  (Lamb  shift),  which  is 
always  present  but,  in  general,  small.  Shifts  and  widths  of  absorption 
lines  are  directly  amenable  to  experimental  observation,  and  hence  a 
verification  of  these  predictions  should  be  feasible. 


IX.  CONCLUSIONS 


Coherences  play  a  pronounced  role  in  the  relaxation  of  adsorbates  if 
the  widths  of  the  transitions  (relaxation  constants)  become  comparable  to  or 
exceed  the  level  separations.  In  the  transient  region  the  time  evolution  of 
the  populations  of  the  adsorbate  states  couples  to  the  evolution  of  the 
coherences,  which  prohibits  a  description  in  terras  of  a  master  equation. 

For  t  -*  o°  (steady  state,  thermal  equilibrium)  the  coherences  vanish,  and  a 
t ine- independent  master  equation  emerges,  which  contains  the  Golden-Rule 
transition  constants  as  parameters.  Measurable  quantities  which  are 
determined  by  a  two-time  quantum  expectation  value,  like  correlation 
functions  or  spectrally- resolved  properties,  however,  involve  the  time- 
evolution  operator  -  if  for  the  density  matrix  p^(t).  This  includes  all 
couplings  between  populations  and  coherences,  and  between  coherences  among 
each  other.  Even  in  the  steady  state,  where  the  coherences  disappear,  their 
tine-regression  operator  obviously  does  not.  Therefore,  a  correct 
evaluation  of  steady-state  properties  requires  that  the  coupling  with 
coherences  is  retained.  In  other  words,  a  random-phase  or  secular 
approximation  is  not  exact  in  general. 

We  have  applied  a  Liouville-operator  formalism  to  derive  a  condensed 
form,  Eq.  (3.3),  of  the  equation  of  motion  for  the  reduced  density  operator 
Pg(t)  of  the  adsorbate.  After  some  algebraic  manipulations,  and  an 
expansion  in  matrix  elements,  we  arrived  at  Eq .  (5.2),  which  generalizes  the 
master  equation  (1.1).  It  should  be  stressed  that  Eq.  (5.2)  contains  the 
sane  parameters  as  Eq.  (1.1),  so  that  no  additional  information  about  the 
system  is  required.  Only  the  coupling  between  the  different  matrix  elements 
is  more  complicated.  Discarding  various  non-secular  couplings  yields  again 
Eq.  (1.1). 


Subsequently  we  have  considered  the  probe-absorption  spectrum  I(uj), 
which  was  expressed  in  the  resolvent  (w  -  +  iT)  ^  in  Eq.  (7.9).  With  a 

two- level  example  it  has  been  illustrated  how  the  formal  expression  can  be 
evaluated,  and  what  the  relevance  of  the  coherence-coherence  coupling  can  be 
for  the  spectral  distribution.  The  results  have  been  compared  with  the 
secular  approximation,  where  I(oj)  is  a  Lorentzian. 
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CAPTIONS 


Fig.  1.  Absorption  profile  l(w)  from  Eq.  (8.9)  as  a  function  of  w.  The 


.-1 


overall  factor  *  ai2^a21  +  a12^  *s  C  en  t0  un^ty- 


Pig.  2. 


Frequencies  are  in  units  of  u)q,  and  the  relaxation  constant  equals 
H  =  0.5.  The  profile  which  is  symmetric  around  is  the  secular 
approximation  I(to)s  from  Eq.  (8.11). 

Same  as  Fig.  1.  but  with  n  *  1.5.  In  this  case  the  width  of  the 
transition  is  larger  than  the  level  separation,  and  it  is  seen  that 


the  resonance  at  vanishes  completely.  The  combined  system  of 


adsorbate  and  substrate  gives  rise  to  a  resonance  near  uj  =  0,  which 
is  not  found  in  the  secular  limit.  There  the  peak  at  u)q  is  smeared 
out  to  form  a  continuous  background.  This  exhibits  clearly  the 
significance  of  the  non-secular  couplings  if  the  levels  have 


overlapping  resonances. 
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